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XMM-JVEWTOJV REVEALS THE QUASAR OUTFLOW IN PG 1115+080 
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ABSTRACT 

We report on an observation of the Broad Absorption Line (BAL) quasar PG 1115+080 performed with 
the XMM-Newton observatory. Spectral analysis reveals the second case of a relativistic X-ray absorbing 
outflow in a BAL quasar. The first case was revealed in a recent observation of APM 08279+5255 with 
the C/ianiira X-ray Observatory. As in the case of APM 08279+5255, the observed flux of PG 1115+080 is 
greatly magnifled by gravitational lensing. The relatively high redshift {z — 1.72) of the quasar places the 
redshifted energies of resonant absorption features in a sensitive portion of the XMM-Newton spectral 
response. The spectrum indicates the presence of complex low-energy absorption in the 0.2-0.6 keV 
observed energy band and high-energy absorption in the 2-5 keV observed energy band. The high- 
energy absorption is best modeled by two Gaussian absorption lines with rest-frame energies of 7.4 keV 
and 9.5 keV. Assuming that these two lines are produced by resonant absorption due to Fe xxv Ka, we 
infer that the X-ray absorbers are outflowing with velocities of ~ 0.10c and ~ 0.34c, respectively. We 
have detected significant variability of the energies and widths of the X-ray BALs in PG 1115+080 and 
APM 08279+5255 over timescales of 19 and 1.8 weeks (proper-time), respectively. The BAL variability 
observed from APM 08279+5255 supports our earlier conclusion that these absorbers are most likely 
launched at relatively small radii of < 10^^(A/B/f /Afg)^/^ cm. A comparison of the ionization properties 
and column densities of the low-energy and high-energy absorbers indicates that these absorbers are 
likely distinct; however, higher spectral resolution is needed to confirm this result. Finally, we comment 
on prospects for constraining the kinematic and ionization properties of these X-ray BALs with the next 
generation of X-ray observatories. 

Subject headings: galaxies: active — quasars: absorption lines — quasars: individual (PG 1115+080) 
— quasars: individual (APM 08279+5255) — X-rays: galaxies — gravitational lensing 



1. INTRODUCTION 

Most of our current understanding of the kinematic 
structure in Broad Absorption Line quasar (BALQSO) 
outflows stems from numerous studies of the velocity struc- 
tures of BAL features that appear blueward of resonant 
UV emission lines. The similarities between the emission- 
line properties of quasars with and without BAL features 
led Weymann et al. (1991) to conclude that BALQSOs 
are "normal" quasars with a BAL global covering factor 
of ~ 10%, viewed along lines of sight that traverse the out- 
flowing UV absorbers. In this model, quasar winds exist 
in most quasars; however, because of the relatively small 
opening angles of these outflows only a small fraction of 
radio-quiet quasars have detectable BAL features in their 
UV and/or optical spectra. Spectropolarimetric observa- 
tions, the statistics of gravitationally lensed BALQSOs, 
improved assessments of incompleteness, and theoretical 
arguments suggest that present optical surveys may be bi- 
ased against the detection of BALQSOs (e.g., Goodrich 
1997; Chartas 2000; Krolik & Voit 1998; Tolea, Krohk, & 
Tsvetanov 2002; Hewett & Fohz 2003). These analyses 
indicate that the true fraction of BALQSOs and BAL cov- 
ering factors may be substantially larger ( > 20%) than 
the often quoted value of ^ 10%. 

The study of BALQSOs in the X-ray band may provide 
valuable information regarding the total column densities 
and the ionization equilibrium of the BAL absorber. X- 



ray absorption in BALQSOs is expected to be largely due 
to photoelectric, bound-free absorption and is less sensi- 
tive to effects such as saturation and the degree of ion- 
ization, often encountered in optical and UV observations 
(e.g., Arav 1998). Because of the heavy X-ray absorption 
in BALQSOs, the X-ray spectra obtained to date have 
modest (e.g., PG 2112+059) to poor signal-to-noise (S/N) 
ratios and cannot accurately constrain the column densi- 
ties, ionization state, geometry, or kinematic properties of 
BAL absorbers (e.g., Gallagher et al. 2001; Green et al. 
2001; Gallagher et al. 2002). Furthermore, it is not clear 
whether the X-ray and UV absorbers are distinct. 

Our recent Chandra observation of the BALQSO 
APM 08279+5255 indicated that quasars may contain 
winds with mass outflow rates considerably larger than 
previously thought (Chartas et al. 2002). This was sug- 
gested by the identiflcation of two X-ray absorption fea- 
tures at rest-frame energies of 8.1 and 9.8 keV. By in- 
terpreting these features as resonant absorption lines of 
highly ionized Fe (e.g., Fe xxv Ka), we inferred the pres- 
ence of X-ray BAL material outflowing from the central 
source at velocities of ~ 0.2c and ~ 0.4c. Our analysis of 
APM 08279+5255 also suggested that the X-ray BAL ma- 
terial in this object is launched at radii of < 2 x 10^^ cm 
and therefore is likely not associated with the UV BAL ma- 
terial which is thought to be launched at larger radii. An 
observation of APM 08279+5255 performed with XMM- 
Newton ~ 1.8 weeks (proper-time) after the Chandra ob- 
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servation showed an absorption feature near 7.4 keV (rest- 
frame) that was interpreted by Hasingcr, Schartel, & Ko- 
mossa (2002, hereafter H02) as a highly ionized iron edge 
(Fe XV-Fe XVlll). If this interpretation is accepted an over- 
abundance of Fc (factor ~ 2 5) is required to fit the XMM- 
Newton spectrum, assuming that the observed low and 
high-energy absorbers are the same. Recently, an ionized 
outflow with a velocity of 0.08-0. Ic has also been detected 
in XMM-Newton observations of the narrow emission-line 
quasar PG 1211+143 (Poimds ct al. 2003). Interestingly, 
there are no detected UV BALs from PG 1211-M43. Ad- 
ditional insight into the association between the X-ray and 
UV BAL material was recently provided by contemporane- 
ous Chandra and HST STIS observations of the BALQSO 
PG 2112-1-059 (Gallagher et al. 2003). Significant differ- 
ences in the variability of the X-ray and UV absorbers in 
PG 2112-1-059 suggest the distinct nature of the X-ray and 
UV BAL material. 

The confirmation of rclativistic X-ray absorbing out- 
flows in most quasars would imply higher mass-outflow 
rates than those based on outflow velocities derived from 
UV BALs. Such outflowing winds could provide a signif- 
icant contribution to the enrichment of the vicinities of 
the quasars and their host galaxies with high-metallicity 
material. To address the question of whether relativistic 
X-ray absorbing outflows are a common phenomenon in 
quasars, we performed a deep observation of PG 1115-1-080 
{z = 1.72, B = 16.1) with XMM-Newton. PG 1115+080 is 
a well-suited object for probing the kinematics of quasar 
winds. The relatively small velocity spread observed for 
the UV absorption in PG 1115+080, compared to the typ- 
ical range of 5,000 - 25,000 km s'^ observed in BAL 
quasars, suggests that PG 1115+080 be classified as a 
mini-BALQSO (Turnshek 1988; Barlow, Hamann, & Sar- 
gent 1997). 

The large flux magnification, estimated to be a factor 
of ~ 20-40 (e.g., Impey et al. 1998), produced by gravi- 
tational lensing provided a bright background source that 
allowed the detection of absorption features from the out- 
flowing wind. Here we present constraints on the kine- 
matic and ionization properties of the outflowing wind of 
PG 1115+080 and provide results from our investigation 
of variabUity of the X-ray BALs of PG 1115+080 and 
APM 08279+5255. Throughout this paper we adopt a 
A-dominated Universe with Hq = 70 km s~^ Mpc~^, f^A 
= 0.7, and Q.m = 0.3. 

2. OBSERVATIONS AND DATA ANALYSIS 

PG 1115+080 was observed for 62.6 ks on 2001 Nov 25 

with XMM-Newton. It was also observed with Chandra 
on 2000 June 2 and 2000 Nov 3 for 26.8 ks and 10.0 ks, 
respectively. The spectral analysis of the Chandra ob- 
servations were presented in Gallagher et al. (2002). In 
§3.2, we provide a comparison between the XMM-Newton 
and Chandra spectral results and include an update on 
the Chandra analysis of PG 1115+080 that incorporates 
a recently developed contamination model for the ACIS 
instrument. 

No significant background flares were present during 
the XMM-Newton observation, and the average back- 
ground rates for the European Photon Imaging Cam- 
era (EPIC) PN and MOS detectors in the 0.2-10 keV 



band were 4.4 x 10~^ events s~^ arcsec"^ and 8.3 x 
10""^ events s~-^ arcsec"^, respectively. We analyzed 
the data with the standard analysis software SAS 5.3. 
PG 1115+080 was clearly detected with the PN, MOSl, 
and M0S2. The total number of source counts detected in 
the 0.2-10 keV band was ~ 12,200 counts for the PN detec- 
tor and ~ 8,185 counts for the combined MOSl and M0S2 
detectors (henceforth referred to as MOSl+2). Flux vari- 
ability was detected at the > 90% confidence level based 
on a Kolmogorov-Smirnov (K-S) test. The timcscalc of 
the flux variations is of order 1000 s. A description of 
the timing analysis and inferred constraints on the Hubble 
constant based on the estimated time delay between the 
lensed images is presented in a separate paper (Chartas et 
al. 2003). We filtered the PN and MOS data by selecting 
events with PATTERNS in the 0-4 and 0-12 ranges, respec- 
tively. The MOSl and M0S2 event lists were combined 
using the SAS task merge to increase the S/N ratio of the 
extracted spectrum. 

The extracted spectra of PG 1115+080 from the PN and 
MOS were grouped to obtain a minimum of 100 counts in 
each energy bin, allowing use of statistics. Background 
spectra for the PN and MOSl+2 detectors were extracted 
from source- free regions near PG 1115+080. The PN and 
MOSl+2 spectra were then fitted simultaneously with a 
variety of models employing XSPEC version 11.2 (Arnaud 
1996). We tested the sensitivity of our results to the se- 
lected background and source cxtrac;tion regions by vary- 
ing the locations of the background regions and varying 
the sizes of the source extraction regions. We did not 
find any significant change in the background-subtracted 
spectra. For all models we included Galactic absorption 
due to neutral gas with a column density of Nii=3.53 x 
10^° cm~2 (Stark et al. 1992). AU quoted errors are at 
the 90% confidence level unless mentioned otherwise with 
all parameters taken to be of interest except absolute nor- 
malization. 

We first simultaneously fitted the PN and MOSl+2 
spectra of PG 1115+080 with a model consisting of a 
power law with neutral intrinsic absorption at z = 1.72 
(fit 1 of Table 1). This fit supports the presence of an 
intrinsic absorber with a column density of Nii= 2.7~to[5 
X 10^^ cm~^. The fit is not acceptable in a statistical 
sense [reduced = 1-54 for 209 degrees of freedom (dof)]. 
The fit residuals show significant absorption from 2-5 keV 
that contributes to the unacceptable fit. To illustrate the 
presence of these features and the absorption features be- 
low 0.6 keV, we flt the spectrum between 1-2.6 keV and 
5.2 10 keV with a power- law model modified by Galac- 
tic absorption and extrapolated this model to the energy 
ranges not fit (see Figure la). For clarity we only show 
the higher S/N ratio PN data in Figure 1; however, all fits 
were performed simultaneously to the PN and MOSl+2 
data. In Figures 2a and 2b we show the Ax residuals be- 
tween the best-fit absorbed power-law model and the PN 
and MOSl+2 spectra. The complex absorption features 
from 2- 5 keV are detected in both the PN and MOSl+2 
spectra. 

We proceeded to model the 2-5 keV features with a vari- 
ety of models. We first considered an absorbed power-law 
model with two cosmologically redshifted Gaussian ab- 
sorption lines near the absorption features appearing at 
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energies of 2.7 and 3.6 keV (observed-frame). This fit is 
a significant improvement compared to the previous one 
at the > 99.99% confidence level and yields a reduced 
= 1.33 for 201 dof (fit 2 of Table 1). The addition of the 
two Gaussian lines significantly reduced the residuals, and 
the remaining most significant contribution to the large 
value of arises from the residuals below 0.6 keV. These 
low-energy residuals are commonly detected in moderate 
S /N ratio X-ray spectra of B ALQSOs and are thought to 
arise from ionized or partially covering intrinsic absorp- 
tion (e.g., Gallagher et al. 2002). To test this we re- 
placed the neutral absorber in our spectral model with 
an ionized intrinsic absorber (see fit 5 of Table 1). In 
particular, we used the absori model contained in XSPEC 
(Done et al. 1992). In Figure lb we show the PN spec- 
trum fit with this model and the residuals of this fit. We 
found a significant improvement in fit quality at the > 
99.99% confidence level (according to the i^-test) with re- 
duced x^ = 1-16 for 200 dof (see fit 5 of Table 1). The 
best-fit rest-frame energies and widths of the absorption 
features are Eabsi = 7.381q }q keV, aabsi < 0.44 keV and 
Eabs2 = 9.5i;i keV (68% confidence), aabs2 = 2Atli keV 
(68% confidence). ^ 

Both of the 2-5 keV features occur in a well-calibrated 
energy region where the effective area of the combined 
telescope and instrument responses varies smoothly. The 
best-fit value for the photon index is F = 1.90to;o4) consis- 
tent with those of typical radio-quiet quasars. The ionized 
absorber is characterized by an ionization parameter of 
^ = L/nr^ = 101.4^43 erg cm s^^ and a hydrogen column 

density of Nu = l.ltoi ^ 10^^ cm~^, where L is the inte- 
grated 5 eV-300 keV incident luminosity, n is the number 
density of the absorber, and r is the distance between the 
absorber and ionizing source. 

We also modeled the absorption from 2-5 keV with one 
and two absorption edges. The fit with one edge is slightly 
worse (at the 90% confidence level) than the fit that in- 
cluded the two Gaussian lines and yields a reduced x^ 
= 1.18 for 206 dof (fit 3 of Table 1). The best-fit edge 
energy and optical depth are Eedge = ^.OTIqJs keV and 
Tedge = 0.31 ± 0.12. An iuspectiou of this fit shows 
that significant residuals are present from 2 5 keV. The 
fit with two absorption edges is comparable in quality to 
the fit that included two Gaussian absorption lines and 
yields a reduced x^ = 1-16 for 204 dof. The best-fit en- 
ergies and optical depths of the two absorption edges are 

Eedgel = 7 Mt^oH keV, Tedgel = 0.31 ±0.11 and Eedge2 

= QMIIoaI keV, Tedge2 = 0.22 ± 0.14. The best-fit val- 
ues for the edge energies would imply that the first edge 
arises from nearly neutral Fe with a 2(t confidence range 
from Fe I to Fe Vii, whereas the second edge energy lies 
Icr above the threshold energy of the highest ionization 
level of Fe xxvi (9.278 keV). While the fit that incorpo- 
rates two absorption edges is acceptable in a statistical 
sense, the inferred energies and optical depths are difficult 
to explain with a physically acceptable model (see §3.1). 
The multiple-edge model would imply that the detected 
high-energy absorption is produced by two distinct ab- 
sorbers, one nearly neutral and the other highly ionized. 

^ The unusually large width of the absorption line at 9.5 keV may 
spectra are needed to resolve this issue. 



The optical depths of the edges imply hydrogen column 
densities of ^ 2.7 x 10^^ cm^^ for the edge at 7.03 keV 
and ~ 6.9 x 10^^ cm~^ for the edge at 9.63 keV (for solar 
abundances and adopting the threshold cross sections of 
Verner & Yakovlev 1995). In §3.1 we investigate whether 
spectral models that incorporate direct and scattered radi- 
ation from the central source can provide acceptable fits to 
the low and high-energy absorption features in the XMM- 
Newton spectrum of PG 1115-^080. 

3. DISCUSSION 

Our XMM-Newton observation of PG 1115+080 has re- 
vealed a BALQSO with two absorption features very sim- 
ilar in rest-frame energy to the ones recently identified in 
Chandra observations of the BALQSO APM 08279+5255. 
In §3.1 we investigate three spectral models to explain 
these high-energy absorption features, provide constraints 
on their kinematic and ionization properties, and exam- 
ine a possible connection between the low and high-energy 
absorbers. In §3.2 we investigate possible variability of 
these absorption features by comparing the Chandra and 
XMM-Newton observations of PG 1115+080 separated by 
19 weeks (proper-time). We also briefly discuss the vari- 
ability of similar absorption features detected in Chandra 
and XMM-Newton observations of APM 08279+5255 sep- 
arated by 1.8 weeks (proper-time) and place constraints 
on the size and location of these absorbers. 

3.1. Origin and Properties of the X-ray Absorbers in 
PG 1115-h080 

We attempted to determine the origin of the low and 
high-energy absorption features seen from PG 1115+080 
by considering the three following models: 

(a) Our first model consisted of neutral and ionized ab- 
sorption with variable metal abundances. The neutral 
component consisted of a power law modified by the XSPEC 
model zvphabs which assumes photoelectric absorption 
with variable abundances. The ionized component con- 
sisted of a power law modified by the XSPEC model absori. 
The Fe abundances of both components were allowed to 
vary during the spectral fits. This multi-component model 
was an attempt to account for the results of the spec- 
tral fits that assumed a two-edge model (described in §2) 
to fit the high-energy absorption features. We assumed 
that our lines of sight through the neutral and ionized ab- 
sorber are different. This can occur if the observed spec- 
trum consists of (1) direct absorbed and (2) scattered, ab- 
sorbed emission from the central source. The fit with this 
multiple-component model was acceptable in a statistical 
sense yielding a reduced x^ = 1-13 for 202 dof. However, 
there are problems with explaining the resulting values for 
several of the parameters. First, the best-fit value of the Fe 
abundance needed to model the neutral absorption feature 
near ~ 7.1 keV is ~ 425j^||Q times solar, and the Fe abun- 
dance needed to model the ionized low-energy absorption 
is ~ 15 ± 10 times solar (errors are at the 90% confidence 
level). These values are unphysically large. Second, the 
Ax residuals of this fit indicate that the double line of 
sight model can fit the low-energy absorption well, but, 
significant residuals remain near the absorption feature at 

due to multiple unresolved narrower components. Higher resolution 
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9.6 keV. We considered the case where the ionized and 
neutral absorbers were along the same line of sight. The 
fit with this model was not acceptable in a statistical sense 
with a reduced = 1-26 for 205 dof. The best-fit value of 
the column density of the neutral absorber for this model 
is A^H = O.Ol^:J X lO^i cm-2. 

(b) We also investigated whether intervening absorbers 
in the lens galaxy at z = 0.31 or in possible damped Ly- 
man alpha systems along the line of sight could explain 
the high-energy absorption features. Spectroscopic obser- 
vations of PG 1115-1-080 in the UV have not detected any 
Lyman-limit discontinuities (Green et al. 1980; Tripp et al. 
1990) and place an upper limit on the lens galaxy's column 
density of neutral hydrogen of A^h J;3 x 10^^ cm^^. The 
most abundant elements with K-shell threshold energies 
near the energies of the observed high-energy absorption 
features are Si and S. For the purpose of this investiga- 
tion we compared the K-shcll threshold energies of Si and 
S (corrected for the redshift of any intervening absorber 
between us and PG 1115-1-080) to the best-fit energies of 
Eedgei = 2.59 keV and Eed,ge2 = 3.54 keV (observed-frame) 
obtained from our fits that assumed two absorption edges 
(see §2 and fit 4 of Table 1). The energies of the Si i 
- XIV and S I - XVI K-shell threshold energies range be- 
tween 1.846-2.673 keV and 2.477-3.494 keV, respectively. 
We expect most of the Si and S in intervening damped 
Lyman alpha systems or in the lens galaxy to be neutral 
or mildly ionized. Wc; find that all the K-shcll threshold 
energies of Si and S fall significantly below the energies of 
the observed high-energy absorption features for absorp- 
tion produced in the lens galaxy with the exception of S 
XVI which falls between the two observed high-energy ab- 
sorption features. All ionization levels of Si and S are also 
ruled out as producing the absorption edge at 3.54 keV for 
any intervening absorber. We conclude that the observed 
high-energy absorption features cannot be produced by 
absorption in intervening systems or in the lens galaxy, 
and the most likely origin is intrinsic absorption by highly 
ionized iron. 

(c) Finally, we examined a model that assumes the high- 
energy absorption features arise from different zones of 
a highly ionized outflowing wind. A similar model was 
used to explain the high-energy absorption features in 
the Chandra spectrum of APM 08279-1-5255. In partic- 
ular, for the case of APM 08279-1-5255 we interpreted 
these features as two distinct absorption systems in the 
quasar outflow with velocities of ~ 0.2c and ~ 0.4c, re- 
spectively; however, see H02 for a different interpreta- 
tion of the XMM-Newton spectrum of APM 08279-^5255. 
For the XMM-Newton observations of PG 1115-^080 we 
found that the absorption features arc located at rest- 
frame energies of 7.38^0' keV (90% confidence) and 
9.50l'j';^o kcV (68% confidence); sec fit 5 of Table 1. Ap- 
plying similar arguments to infer the nature of the lines in 
PG 1115-h080 as those presented in Chartas et al. (2002) 
for APM 08279-1-5255, we conclude that the absorption 
lines detected from PG 1115+080 are most likely due to 
highly ionized Fe ions in the quasar wind. We estimate 
that the 7.38 keV and 9.50 keV absorption features cor- 
respond to wind velocities (depending on the ionization 
state of Fe) of 0.10c (Fe xxv Ka), 0.04c (Fe xxvi Ka) 
and 0.34c (Fe xxv Ka), 0.30c (Fe xxvi Ka) respectively, 



relative to the systemic redshift of z = 1.72. For these 
velocity calculations we included the special relativistic 
velocity correction. 

At the la level, the ratio of the two absorption-line en- 
ergies {Eabsi/ Eabs2) IS inconsistcut with the ratio of the 
energies of the iron Ka to K/3 transitions of Fe xxv and 
Fe XXVI. The predicted ratio of EWabsi/ EWabsi of the 
two high-energy absorption lines is ~ 7.1 (7.4) assum- 
ing that these lines are the iron Ka and K/3 transitions 
of Fe XXV (Fe xxvi) and the absorption is on the lin- 
ear part of the curve of growth. The observed ratio of 
EWabs\/EWabs2 = 0.1 ± 0.06 is much less than the pre- 
dicted ratio, and thus docs not obviously support the iron 
Ka and K/3 interpretation of the high-energy absorption 
lines. We note, however, that in the likely case where the 
absorption lines are saturated or are composed of multi- 
ple unresolved components the observed EWabs2 / EWabsi 
ratio could deviate significantly from the predicted ratio 
assuming the iron Ka and K/3 interpretation. A second 
argument against the iron Ka and K/3 interpretation is 
the observed radical difference of the line profiles of the 
absorption lines at 7.38 keV and 9.50 keV (see fit 5 of 
Table 1) that is not expected if these absorption features 
are produced by iron Ka to K(3 transitions of Fe xxv or 
Fe xxvi. 

UV spectra of PG 1115-1-080 acquired with the Interna- 
tional Ultraviolet Explorer {WE) reveal an O vi A1033 
resonant emission line accompanied by absorption in the 
line core and absorption features blueward of the line that 
extend over a velocity range of 0-6000 km s~^ (Michalit- 
sianos & Oliverson 1996). A possible explanation for the 
large differences between the UV and X-ray BAL wind 
velocities is that the two absorbers are launched at dif- 
ferent distances from the central ionizing source and are 
therefore accelerated to different terminal velocities. To 
examine this hypothesis further wc; estimated the wind 
velocity as a function of radius following §5 of Chartas et 
al. (2002). We assumed a minimum launching radius for 
the UV BAL absorber of r„un 1 x \Q^^{Mbh/Mii)^l'^ cm, 
based on hydrodynamical simulations performed by Proga 
et al. (2000), where Mbh is the mass of the black hole of 
PG 1115-h080 and Mg = IO^Mq. In Figure 3 we plot the 
velocity of the UV absorber versus radius from the central 
source for launching radii of rmin, 2.5rniin and 5rinin- We 
have assumed a force multiplier of F/ = 100 (e.g., Arav, 
Li, & Begelman 1994), Luv = 5 x lO**** erg s~^, LboI = 
2 X 10"^^ erg s~^, and Lsoi/LEdd = 0.1. The terminal ve- 
locity of UV material launched at the minimum launching 
radius is ~ 0.07c, significantly less than the velocities of 
the X-ray BAL absorbers. As wc proposed in the case of 
APM 08279-f5255, a possible explanation for the larger 
velocities of the X-ray BAL absorbers is that this material 
is initially accelerated at smaller launching radii than the 
UV BAL material. The short recombination timescales of 
Fe XXV and/or Fe xxvi compared to the estimated time 
for the X-ray BAL material to be accelerated close to its 
terminal velocity imply that the observed highly ionized 
Fe is located near the launching radius where it would still 
be exposed to intense ionizing radiation from the central 
X-ray source. One of the implications of small launch- 
ing radii for the X-ray BAL material is possible variability 
over timescales of weeks to months of the energies and 
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widths of the X-ray BALs as the ionized X-ray absorbers 
are accelerated to their terminal velocities. 

We investigated a possible connection between the low 
and high-energy X-ray absorbers by comparing their kine- 
matic, ionization, and absorbing properties. We searched 
for a possible outflow velocity of the low-energy absorber 
by allowing the redshift parameter of the ionized-absorber 
model to be free. The best-fit redshift of the absorber 
of z = 1.53^0 7 (90% confidence errors) was consistent 
with the systemic value of z — 1.72 (and also with the 
blue-shifted values of z = 1.46 and 0.92 expected from 
the velocities of the high-energy absorbers), and we found 
no significant improvement in the low-energy fit residuals 
and no improvement in the global fit. We note, however, 
that if the low-energy absorber has a large velocity dis- 
persion it would be difficult to detect a blue-shift. The 
estimated hydrogen column density of the low-energy ab- 
sorber ranges from A^h=1-8 x 10^"'^ cm~^ assuming the 
neutral absorption model (fit 2 of Table I) to A^h^I-I x 
10^^ cm^^ assuming the ionized absorber model (fit 5 of 
Table 1); we note that these A^h values may have signif- 
icant systematic errors due to the inadequacy of the ab- 
sorption models. Specifically, these models do not include 
effects such as possible velocity dispersion of the absorber, 
BALs, and multiple ionization zones. To obtain an esti- 
mate of the hydrogen column densities of the high-energy 
absorbers we performed a simple curve-of-growth analysis 
(Spitzcr 1978) on the detected absorption features. For 
no ionization correction (assuming all of the Fe is in the 
Fe XXV ionization state) and assuming solar abundances, 
the total hydrogen column densities corresponding to the 
absorption lines at 7.38to\[l keV and 9.50t°:4S keV are 
iVn = 4 X 10^^ cm^^ and A'h = 4 x 10^'^ cm^^, respec- 
tively. We consider these estimates of the column densities 
also to have likely systematic errors since the absorption 
lines may contain multiple unresolved absorption compo- 
nents leading to an overestimate of the column densities 
of individual velocity components and the b parameters 
obtained from fits of Gaussian lines to observed absorp- 
tion features are valid only in the optically thin limit. 
On the other hand, possible saturation of the absorption 
lines will lead to an underestimate of the column densities. 
Our preliminary comparison of the kinematic properties of 
the low-energy and high-energy absorbers cannot conclu- 
sively show whether they are outflowing at similar veloc- 
ities, however, our spectral analysis indicates that these 
absorbers have different ionization parameters. In partic- 
ular, the best-fit ionization parameter for the low-energy 
absorption is ^ = 101^^48 erg cm s^-'^ (sec fit 5 of Table 1), 
whereas, the presence of highly ionized Fe (high-energy 
absorption features) requires ^ > 1000 erg cm s^^. We 
therefore conclude that the low and high-energy absorbers 
of PG 1115+080 are likely distinct; however, higher spec- 
tral resolution is needed to confirm this result. 

3.2. Variability of X-ray BALs in PG 1115+080 and 
APM 08279+5255 

We searched for variability of the X-ray BALs of 
PG 1115+080 and APM 08279+5255 by comparing their 
properties at different epochs. We first compared the X- 

* ACISABS is an XSPEC model contributed to the Chandra users so 
soft/ soft-list . cgi 



ray spectra of PG 1115+080 obtained in Chandra and 
X MM- Newton observations separated by ~ 19 weeks. We 
fit the Chandra spectrum obtained on 2000 June 2 for 
26.8 ks simultaneously with the XMM-Newton spectrum 
of PG 1115+080 using the spectral model of fit 5 of Ta- 
ble 1. To account for the recently observed quantum effi- 
ciency decay of ACIS, possibly caused by molecular con- 
tamination of the ACIS filters, we have applied a time- 
dependent correction to the ACIS quantum efficiency im- 
plemented in the XSPEC model ACISABS 1 . 1. We allowed 
the normalization of the power-law component to be a free 
parameter in the fit to account for any flux variability 
and/or systematic calibration uncertainties between the 
Chandra and XMM-Newton instrumental efficiencies. We 
found a degradation in fit quality at the 90% confidence 
level compared to the fit of the same model to the XMM- 
Newton spectrum of PG 1115+080 based on the F-test. 
This change in fit quality is suggestive of variability of the 
spectrum of PG 1115+080. 

When we repeated the fit to the Chandra spectrum of 
PG 1115+080 allowing the model parameters of the low 
and high-energy absorbers to vary, we obtain best-fit rest- 
frame line energies, widths, equivalent widths, and 68% 
confidence levels of the absorption features of Eabsi = 
7.37i°i^ keV, aabsi < 0.35 keV, EWabsi = eV 
and Eabs2 = 8.5toi keV, aabs2 < 0.25 keV, EWabs2 = 
178^92^ eV. The best-fit values for the ionization parame- 
ter and hydrogen column density of the low-energy ionized 
absorber during the Chandra observation are ^ = 398^2io 
erg cm s~^ and A^h = 3.2l^;5 x 10^^ cm~^. A com- 
parison between the Chandra and XMM-Newton obser- 
vations indicates that the line width of the second absorp- 
tion feature varied to <Jabs2 — 2.4^q' i keV for the XMM- 
Newton observation. The properties of the absorption line 
at Eabsi = 7.4 keV were consistent between both obser- 
vations. In Figures 4a and 4b we show the Ax residu- 
als between the best-fit simple absorbed power-law model 
and the XMM-Newton PN and the Chandra ACIS spectra 
of PG 1115+080, respectively. The ionization parameter 
and hydrogen column density of the low-energy ionized ab- 
sorber also appeared to vary to ^ = 101 tig erg cm s~^ and 
A^H = l-ltoi X 10^2 cm-2 for the XMM-Newton obser- 
vations. In Figures 5a and 5b we show the 68% and 90% 
confidence contours of photon index versus intrinsic low 
energy absorption column density and ionization parame- 
ter versus intrinsic column density for the XMM-Newton 
and Chandra observations assuming fit 5 of Table 1. These 
contours indicate a possible change (at the 68% confldence 
level) of the intrinsic column density and the ionization 
parameter between observations. To illustrate the signifi- 
cance of variability in the high-energy absorption features 
detected in PG 1115+080 with XMM-Newton and Chan- 
dra wc present in Figures 6a and 6b confidence contours 
of £",15.52 vs. (Tabs2 and Eabs2 vs. line-fiux, respectively. 
We note, howcvcir. that these confidence contours are only 
suggestive of variability in the width and line flux of the 
absorption line at Eabs2 since this line is only detected 
at the 68% confidence level in the Chandra spectrum of 
PG 1115+080. Flux and possibly spectral variability in 

3 exchange world wide web-site http://asc.harvaj:d.edu/cgi-gen/cont- 
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the 0.2 - 1.0 kcV band have boon detected in earlier obser- 
vations of this object performed with the Einstein IPC and 
ROSAT PSPC on 1979 Dec 10 and 1991 Nov 21, respec- 
tively (see Figures 5 and 6 of Chartas 2000). Specifically, 
the 0.2 2 kcV flux varied from 6.5 x lO^^'' erg s^^ cm~^ 
to 0.4 X 10~^^ erg s~^ cm~^ between the Einstein IPC 
and ROSAT PSPC observations. The 0.2-2 keV and 0.2- 
10 keV fluxes of PC 1115+080 observed with the PN are 

2.7 X 10~^^ erg s~^ cm~^ and 6.1 x 10~^^ erg s~^ cm~^, 
respectively. 

We also searched for variability of the X-ray absorption 
lines recently detected in Chandra and XMM-Newton ob- 
servations of the BALQSO APM 08279-1-5255. The Chan- 
dra and XMM-Newton spectra were fit with model 5 of 
Table 1. Both fits were acceptable with rc;duc;cd = 1-04 
for 102 dof in the Chandra case and reduced = 1-04 for 
116 dof in the XMM-Newton case. A comparison between 
the two observations separated by ^ 1.8 weeks (proper- 
time) indicates that the values (quoted with 68% confi- 
dence errors) of Eabsi = 8.05to;o4 keV , aabsi < 0.1 keV 
, Eabs2 = ^-81"'! keV, and (Tabs2 = 0.44^;^;} kcV observed 
during the Chandra observations of APM 08279-f-5255 
varied to Eabsi = 8.15tgj| keV, aabsi = 0.6^^;^ keV, 
Eabs2 = ll.Olar keV, and ct„,6,,2 = 2.0^^;^ kcV for the 
XMM-Newton observations. In Figures 7a and 7b we show 
the Ax residuals between the best-fit simple absorbed 
power-law model and the XMM-Newton PN and Chandra 
ACTS spectra, respectively. The detection of significant 
variability of the widths and energies of the high-energy 
absorption lines in APM 08279-1-5255 confirms the intrin- 
sic origin of the high-energy absorption features in this 
BALQSO (e.g., Hamann et al. 1995). The fact that the 
two absorption lines were narrower during the Chandra 
observation than during the XMM-Newton observation of 
APM 08279-1-5255 may explain why these features were 
more easily detectable in the former case. We conclude 
that the observed variability of the energies and widths 
of the X-ray BALs in APM 08279-1-5255 over a timescale 
of 1.8 weeks (proper-time) supports our earlier conclusion 
that these absorbers are most likely launched at relatively 
small radii of less than lO^^iMbh/Ms)^^'^ cm. 

To summarize, we have reported the detection of a rel- 
ativistic outfiow in the BALQSO PC 1115+080 imprinted 
as resonant X-ray absorption lines at rest-frame energies 
oi7.38to.\l keV (90% confidence) and 9.501^0 keV (68% 
confidence). These energies suggest the presence of two 
distinct absorbers outflowing at velocities of ^ 0.10c and 
~ 0.34c. We found significant variability of the properties 
of the X-ray BALs in APM 08279+5255 over a period of 

1.8 weeks (proper-time) supporting our earlier conclusion 
that X-ray BALs are launched at relatively small radii 
compared to the expected UV BAL minimum launching 
radii. 

Resolving the low and high-energy absorption features 

seen from PC 1115+080 and APM 08279+5255 is re- 
quired to constrain better the kinematic, ionization, and 
absorbing properties of the observed outfiows. The X-ray 
calorimeter onboard the ASTR0-E2 satellite (scheduled to 
launch in early 2005) has an effective area that is » 4 times 
larger than that of the Chandra gratings in the energy 
range of the observed X-ray BALs and an energy resolu- 
tion of about 12 eV. Future observations of PC 1115+080 



and APM 08279+5255 with ASTR0-E2 may provide high- 
resolution spectra of X-ray BALs with exposure times of » 
300 ks. The current plan for the Constellation — X obser- 
vatory calls for a combination of four satellites that will 
provide the higher energy resolution and larger effective 
area that is needed to resolve the complex structure of the 
X-ray BALs in PC 1115+080 and APM 08279+5255 with 
eflflcient exposure times of about 10 ks. In addition, moni- 
toring of the X-ray BALs with Constellation — X on short 
timescales of the order of days may allow us to track the 
acceleration phase of the absorbers via the velocity shifts 
of the resonant absorption lines. 

We thank those involved with the XMM-Newton Guest 
Observer Facility for assisting in the initial data reduction 
of PC 1115+080. We acknowledge financial support from 
NASA grant NAG5-9949. WNB acknowledges financial 
support from NASA LTSA grant NAG5-13035 and NASA 
grant NAG5-9932. 
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TABLE 1 

RESULTS FROM FITS TO THE XMM-Newton SPECTRUM OF PG 1115+080 



Fit" Model 

1 PL, neutral absorption 
at source. 



Parameter 

r 

JVh 



(0.27 



Value*^ 

1 9i+»"3 

-^■^-^-0.03 

+S:o5) X 10^= 

322.7/209 
7.1 X 10"'' 



PL, neutral absorption 
at source, and two 
Gaussian absorption lines 
at source. 



r 

EWa6,2 



1.76! 



.^ + 0.04 
^-0.04 

(O-lslo.Ss) X 10^2 cm- = 

7A2t°o% keV 
< 0.55 keV (68% errors) 
0.15l°;^? keV 
9.9t°;^ keV 
S.Olog keV (68% errors) 
■ 2.2+_l-l keV 
267.4/201 
1.2 X 10"'* 



PL, ionized absorption 
at source and absorption 
edge. 



r 

JVh 

Tedge 

xV^- 



1 94+o°5 
1.2t°:^) X 10== cr 
103^45 erg cm s~ 



7.07 



+0.20 



-0.23 
+0.12 
-0.11 



keV 



o.si: 

243.2/206 
0.04 



PL, ionized absorption 

at source and two absorption 

edges. 



r 

JVh 

'^e dgel 

1~edge2 

xV^' 
i'CxV^')'' 



, + 0.05 
^-0.02 



1.88T 

(l.otoi) X 10"" cm- 

ni +R0 - 1 

91-46 <=''g ^ 
7.03tn\l kcV 



□ + 0.19 
'-0.32 

'^■■^J^-0.11 



9.63 



+ 0.77 



-0.43 
,+0.14 
-0.14 



keV 



0.22:^ 
236.1/204 
0.06 



PL, ionized absorption 
at source, and two 
Gaussian absorption lines 
at source. 



r 

JVh 

Eaijsl 
faiisl 
EWabal 

Ea6,2 

EW„6=2 

Pix^l^Y 



1 90+°°* 

X 10== cm- 
101+4! erg cm s~^ 
7-38l° iJ keV 
< 0.44 keV 



0.14 



+ 0.06 



-0.06 

9.50l5:4 keV (68% errors) 
2.4+J = keV (68% errors) 
■ 1.4i° s keV 
230.9/200 
0.07 



koV 



"All model fits include Galactic absorption toward the source (Stark et 
al. 1992). All fits are performed on the combined spectrum of images 
Al, A2, B and C of PG 1115+080. 

''All absorption-line parameters are calculated for the rest frame. 
^All errors are for 90% confidence unless mentioned otherwise with all 
parameters taken to be of interest except absolute normalization. 
'^P(x=/i/) is the probability of exceeding %= for v degrees of freedom if 
the model is correct. 
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Fig. 1. — (a) The top panel shows the XMM-Newton PN spectrum of the combined images of PG 1115+080 fit with Galactic 
absorption and a power-law model to events with energies lying within the ranges of 1-2.6 keV and 5.2-10 keV. The lower panel 
shows the residuals of the fit in units of la deviations. Two absorption features from 2.6-5.2 keV are noticeable in the residuals 
plot, (b) The top panel shows the PN spectrum fit with Galactic absorption, ionized absorption at the source, a power-law model, 
and two Gaussian absorption lines. In the lower panel the residuals indicate that this model can account for most of the spectral 
features of PG 1115+080. 



10 



CHARTAS ET AL. 



2.72 



Rest-Frame Energy (keV) 

5.44 13.60 



27.20 



CM - 



CM 
I 



I 



CO 

I 



CM 



CM 
I 



I 



CO 

I 



T 1 1 I 



T 



(a) 



^absl^absS 



MOS 1+2 XMM-Newton 
—I 1 1 1 — I — I— 




(b) 



^absl^abs2 



PN XMM-Newton 



2 5 
Observed— Frame Energy (keV) 



10 



Fig. 2. — Ax residuals between the best-fit Galactic absorption and power-law model and (a) the MOSl+2 spectrum of 
PG 1115-1-080 and (b) the PN spectrum PG 1115-1-080. This model is fit to events with energies lying within the ranges 1-2.6 keV 
and 5.2-10 keV. Eabsi and Eabs2 indicate the best-fit energies of the Gaussian absorption lines obtained in fit 5 of Table 1. 
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Fig. 3. — Wind velocity as a function of radius from the central source of PG 1115+080 for a radiation-pressure driven wind. For 



a qualitative comparison we have estimated the wind velocities for launching radii of 2 x 10^" cm, 5 x 10^" cm, and 1 x lO^'^ cm. 
We have over-plotted the observed O VI BAL (dashed line) and possible Fe xxv BAL (dotted lines) velocities. 
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Fig. 4. — Ax residuals between the best-fit Galaetic absorption and power-law model and (a) the XMM-Newton PN spectrum of 
PG 1115-1-080 and (b) the Chandra ACIS spectrum of PG 1115-1-080. The model has been fit to events with energies lying within 
the ranges of 1-2.6 keV and 5.2-10 keV. Eabsi and Eabs2 indicate the best-fit energies of the Gaussian absorption lines obtained 
using the spectral model described in fit 5 of Table 1. 
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Fig. 5. — 68% and 90% confidence contours between (a) the spectral slope and the intrinsic absorption column density, and (b) 
the ionization parameter and the intrinsic absorption column density for the XMM-Newton (solid contours) and Chandra (dotted 
contours) observations of PG 1115+080 assuming fit 5 of Table 1. 
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Rest-Frame Energy (keV) 



Fk;. 6. (a) 68% and 90% confidence contours between the width of the absorption feature at Eabs2 and the energy Eabs2 for 
the XMM-Newton (sohd contours) and Chandra (dotted contours) observations of PG 1115+080 assuming fit 5 of Table 1. (b) 
68% and 90% confidence contours between the flux of the absorption feature at Eabs2 and the energy Eabs2 for the XMM-Newton 
(solid contours) and Chandra (dotted contours) observations of PG 1115+080 assuming fit 5 of Table 1. These confidence contours 
are only suggestive of variability of the absorption line Eabs2 between the XMM-Newton and Chandra observations since this 
absorption line is only detected at the 68% confidence level in the Chandra spectrum. 
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Fig. 7. — Ax residuals between the best-fit Galactic absorption and power-law model and (a) the XMM-Newton PN spectrum and 
(b) the Chandra ACIS spectrum of APM 08279-1-5255. This model is fit to events with energies lying within the range 2.2-10 keV. 
Eabai and Eabs2 indicate the best-fit energies of the Gaussian absorption lines obtained using the spectral model described in fit 
5 of Table 1. 



